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I. INTRODUCTION 

Rare earth nickelates with the chemical formula 
RMO3, where R stands for a rare earth ion, is one of the 
few families of oxides that undergoes a first order metal- 
insulator phase transition. A metal-insulator transition 
(M-I transition) is an electronic phase transition which 
is usually solid to solid. The first order electronic phase 
transition in many systems is associated with a hystere- 
sis, and they exhibit a phase separated (PS) state in the 
vicinity of the transitional^. The PS state has drawn 
a lot of attention in the case of manganites, mainly due 
to its implications for colossal magnetoresistance, and it 
exhibits a variety of time dependent effects such as relax- 
ation of resistivity^ and magnetization^ and giant 1/f 
noised. 

RNiC>3 crystallizes in the orthorhombically distorted 
perovskite structure with space group P& m „. The de- 
gree of distortion increases with decrease in size of rare 
earth ion as one moves towards the right in the peri- 
odic tabled. The ground state of nickelates (R 7^ La) is 
insulating and antiferromagnetic. On increasing the tem- 
perature these compounds undergo a temperature driven 
antiferromagnetic to paramagnetic transition, and an in- 
sulator to metal transition. For NdNiOa and PrNi03, 
the M-I transition temperature (Tmi) and the magnetic 
ordering temperature (T/v) are the same, while Tmi is 
greater than Tjv for later members of the series such as 
EuNiOs and SmNi03. The Tmi of nickelates increases 
with decreasing size of the rare earth ion and the later 
members of the series also show a charge ordering transi- 
tion associated with the M-I transitional^. The charge 
ordered state has also been observed in epitaxially grown 
thin films of NdNiO^. 

The M-I transition in NdNiOs is associated with a la- 
tent heat and a sudden increase in unit cell volume which 



are signatures of a first order phase transitio n 11 ' 12 . There 
is also a hysteresis in electronic transport properties just 
below Tmt^- The hysteresis in transport measurements 
points toward the coexistence of metallic and insulating 
phases below Tmi^- In this paper we report a detailed 
study of the nature of NdNiOs below Tmi through tem- 
perature and time dependent electrical resistivity mea- 
surements. 



II. EXPERIMENTAL DETAILS 

Polycrystalline NdNiOa samples in the form of 6 mm 
diameter and 1 mm thick pellets were prepared and char- 
acterized as described elsewhere*— The preparation 
method uses a high temperature of 1000°C and a high 
oxygen pressure of 200 bar. 

All the temperature and time dependent measurements 
were done in a home made cryostat. To avoid ther- 
mal gradients in the sample during measurement it was 
mounted inside a thick-walled copper enclosure so that 
during the measurement the sample temperature would 
be uniform. It was found that mounting the sample like 
this improved the reproducibility of the data, especially 
in the time dependence measurements, significantly. A 
Lakeshore Cryotronics temperature controller model 340 
was used to control the temperature and the tempera- 
ture stability was found to be better than 3mK during 
constant temperature measurements. 

Below Tmi NdNiOs is not in thermodynamic equilib- 
rium and slowly relaxes. Because of this the data we get 
depend on the procedure used for the measurement. We 
used the following procedure to measure the temperature 
dependence of resistivity. While cooling we start from 
300 K, and then record the data in steps of 1 K interval 
after allowing the temperature to stabilize at each point. 
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In between two temperature points the sample was cooled 
at a fixed cooling rate of 2 K/min. After the cooling run 
is over we wait for one hour at 82 K and then the heating 
data was collected at every one degree interval. The in- 
termediate heating rate between temperature points was 
the same as the cooling rate used earlier. This cycle of 
measurements was repeated with a different cooling and 
heating rate of 0.2 K/min also. 

It was observed that the resistivity above Tmi does 
not show any time dependence, and it is also indepen- 
dent of measurement history. Thus to avoid the effect 
of any previous measurements, all time dependent ex- 
periments in the cooling run were done as follows: first 
take the sample to 220 K (above Tmi), wait for half an 
hour, then cool at 2.0 K/min to the temperature of inter- 
est and once the temperature has stabilized record the 
resistance as a function of time. In the heating run the 
time dependent resistivity was done in a similar fashion: 
first take the sample to 220 K, wait for half an hour, then 
cool at 2.0 K/min to 85 K, wait for one hour, and then 
heat at 2.0 K/min to the temperature of interest and once 
the temperature has stabilized record the resistance as a 
function of time. 

The four probe van der Pauw method was used to 
measure the resistivity and standard precautions, such 
as current reversal to take care of stray emfs, were taken 
during the measurement. We also took care to ensure 
that the measuring current was not heating up the sam- 
ple. A Keithley current source model 224 and a DMM 
model 196 were used for the resistivity measurements. 

III. RESULTS 

Figure Q] shows the electrical resistivity of NdNiC>3 as 
a function of temperature. We see that the function is 
multiple valued, the cooling and heating data differing 
significantly from each other and forming a large hystere- 
sis loop. The resistivity plot indicates that NdNiC>3 un- 
dergoes a relatively sharp M-I transition at about 190 K 
while heating with a width of about 10 K. In contrast, 
while cooling, the resistivity shows a rather broad M-I 
transition centered around 140 K with a spread of about 
40 K. Below 115 K or so, the heating and cooling data 
merge and the \ogp vs 1/T plot is linear. This indicates 
that the sample is insulating at low temperatures and, 
if the band gap is A, the resistivity should follow the 
relation 

p(T) = p Q exp{A/k B T) (1) 

Below 115 K the data fit quite well to this model, with 
a coefficient of determination, R 2 = 0.99955. po and A 
for the insulating region turn out to be 99 mil cm and 
42meV respectively, which are in reasonable agreement 
with the values previously reported^. 

We collected more hysteresis data with different min- 
imum temperatures such as 140 K, 146 K and 160 K. In 
these measurements we cool the sample from 220 K to 



one of the minimum temperatures mentioned above and 
then heat it back to 220 K, both operations being car- 
ried out at a fixed rate of 2 K/min. Loops formed in this 
fashion are called minor loops and these are indicated by 
the labels a, b and c in Figure [TJ In the cooling cycle 
all the three minor loops coincide with the cooling curve 
of the full hysteresis loop. In the heating cycle, for loops 
a and 6 with lower minimum temperatures, the resistiv- 
ity decreases with increasing temperature and joins the 
full loop at 195 K. In the case of loop c, as we increase 
the temperature, the resistivity increases somewhat till 
about 182 K and then it falls and joins the full loop at 
195 K. 

The resistivity also shows a noticeable dependence on 
the rate of temperature change in the cooling cycle as 
shown in the inset of Figure [TJ The data for the lowest 
curve was collected at 2 K /min and for the middle curve 
at 0.2 K/min. The uppermost curve is an estimate ob- 
tained by extrapolating the time dependence data shown 
in Figure [2] to infinite time. We did not see any rate de- 
pendence in the heating cycle. This is an indication that 
while cooling, below the M-I transition temperature, the 
system is not in equilibrium and hence the resistivity 
evolves with time. The non-rate dependence observed in 
the data while heating tells us that in the heating cy- 
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Figure 1: (Color online) p vs. 1/T plot for NdNi0 3 . The 
blue circles represent cooling data and the red squares stand 
for heating data (cooling/heating rate 2.0 K/min). The solid 
line is a least square fit to the band gap insulator model below 
115 K. The curves labeled a, b & c were taken as described in 
the text. The accuracy of the data points is better than 1% 
everywhere. The error does not exceed 0.8 mficm anywhere. 
The inset shows p vs 1/T for three different cooling rates: 
lower curve (blue stars): 2 K/min, middle curve (red circles): 
0.2 K/min, upper curve (black pluses): infinitely slowly (ex- 
plained in text). The connecting lines are to guide the eyes. 
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Figure 2: (Color online) Time dependence of resistivity while 
cooling, at various temperatures in the range 140 to 155 K, for 
a period of one hour. Maximum time dependence is seen at 
145 K which is about 300%. The curve at the bottom, which 
looks like a straight line, shows the increase in resistivity in a 
heating run taken at 145 K and the change in this case is less 
than 0.2 %. Not all the data are shown here to avoid clutter. 



cle the system is either in, or very close to, equilibrium. 
These observations are corroborated by the data shown 
in Figure El 

A subset of the time dependent resistivity data taken 
while cooling is shown in Figure [21 The data are pre- 
sented as p(t)/p(t = 0) so that the values are normalized 
to unity at t = for easy comparison. We found that 
that below 160 K, the resistivity of the sample increases 
with time considerably. A maximum relative increase 
in resistivity of about 300 % for a duration of one hour 
is seen at 145 K, the time dependence being lower both 
above and below this temperature. We fitted the p(T, t) 
curves in figure 2 to the stretched exponential function 



p(t)=po+pi(l-exp(-(i/r)T)) 



(2) 



where pa, pi, t and 7 are fit parameters. The fits are 
quite good with the R 2 value greater than 0.999 in most 
cases. See Table HI We note that the exponent 7 lies 
in the range 0.5 < 7 < 0.6 and t has a peak around 
147.5 K. The variation of resistivity with time shows that 
the system slowly evolves towards an insulating state at a 
constant temperature. We collected data up to 12 hours 
(not shown here) to check whether the system reaches 
an equilibrium state, but found that it was continuing to 
relax even after such a long time. 

Figure [3] compares the resistivity in a cooling run with 
and without intermediate aging. These data were taken 
as follows: we start from 220 K, come down to 160 K 
at 2K/min, collect time dependence data for one hour 
and after that resume cooling at 2 K/ min and go down 
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Table I: Fit parameters for the time dependence data shown 
in Figure [2] The degrees of freedom of the fits DOF w 1000. 
The x 2 /DOF for 150 K and 155 K are too small, indicating 
that we have overestimated the error in resistivity in these 
cases. Anyway, we note that, the R 2 values are consistently 
good and indicate reasonably good fits. 




125 



135 



T(K) 



145 



155 



Figure 3: (Color online) Temperature dependence of p on 
cooling with ageing by intermediate stops of one hour each 
(blue open circles) at 155, 150, 145, 140, 135, and 130 K. 
Red solid circles show the resistivity without ageing. The in- 
set compares the temperature dependence of pi/ po calculated 
from the fit of Equation (J2]) to p(T, t) data with intermediate 
ageing (blue circles) to the time dependence data shown in 
Figure EJred squares). 



to 155 K, collect resistivity time dependence again for 
one hour and so on down to 110 K for every 5K inter- 
val. We note that when cooling is resumed after ag- 
ing for an hour, p(T) curve merges smoothly with the 
curve obtained without aging. A rather similar observa- 
tion has been reported in the phase-separated manganite 
Lao^Cao.sMno.gsFeo.osO^. 

The red squares in the inset of Figure [3] displays pi/ po 
obtained from stretched exponential fitting of time de- 
pendence resistivity data of Figure[21 The p\ / po increases 
on decreasing the temperature below 160 K, attains a 
maximum around 145 K, and then decreases on further 
lowering the temperature. We have already seen this be- 
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havior in Table |H The value of pi / po obtained from the 
fitting of time dependence data with intermediate aging 
is shown as blue circles in the inset of figure O and we 
note that the value of p\ / po is relatively small compared 
to what one obtains without intermediate aging. 

In the heating cycle, the p vs T exhibits a rather 
sharp insulator to metal transition. The magnitude of 
time dependence in heating runs (maximum w 0.2%) 
is negligible compared to what one gets in cooling runs 
(maximums 300%) (figure [2]), which suggest that in the 
heating run, below Tmi, the sample is almost fully insu- 
lating and stable. 

The major observations that have been made so far are 
summarized below. 

1. The metal to insulator transition in cooling runs 
is rather broad and extends from about 160 K to 
120 K. 

2. In cooling runs a large amount of time dependence 
in resistivity is seen compared to which the time 
dependence seen in heating runs is negligible. 

3. NdNiC>3 undergoes a rather sharp insulator to 
metal transition between 185K and 195 K while 
heating. 

4. There is a large hysteresis between the cooling and 
heating data below the MI transition temperature. 



IV. DISCUSSION 

It has been known that the M-I transition in NdNiC>3 
is associated with a latent heat and a jump in the unit 
cell volume which are the characteristics of a first order 
phase transitio n 11 ' 12 . DSC measurements of Granados et 
al. show that, while cooling, the latent heat released dur- 
ing the M-I transition extends over a broad temperature 
rangeii. This is consistent with the results of the resis- 
tivity measurements that the M-I transition in NdNiC>3 is 
broadened while cooling . Thus we can say that NdNiC-3 
has a broadened first order metal to insulator phase tran- 
sition while cooling. 

We have seen that, below Tmi, while cooling, the sys- 
tem is not in equilibrium and evolves with time which 
suggests that it is in a metastable state. The resistivity 
of this metastable state slowly increases with time which 
means that a slow metal to insulator transition is going 
on in the system. On the other hand while heating up 
from low temperature we saw that the system remains in- 
sulating all the way up to 185 K and then it transitions to 
the metallic state by 195 K. The resistivity was found to 
have negligible time dependence during the whole of the 
heating run. This result indicates that the low tempera- 
ture insulating state is stable and is probably the ground 
state of the system. It is well known that below a phase 
transition temperature a high temperature phase can sur- 
vive as a metastable supercooled state. Based on these 



facts we propose that the metastable state in the case of 
NdNiC-3 consists of supercooled (SC) metallic phases and 
stable insulating phases. 

In a first order transition, a metastable SC phase can 
survive below the first order transition temperature (Tc), 
till a certain temperature called the limit of metastability 
(T*) is reached ^ 14 i 15 . In the temperature range T* < 
T < T c there is an energy barrier separating the SC 
phase from the stable ground state. The height of the 
energy barrier (U) separating the SC metastable phase 
from the stable ground state can be written as U oc f(T— 
T*), where / is a continuous function of (T — T*), and 
vanishes for T < T*. As the temperature is lowered, at 
T = T*, the SC metastable phase becomes unstable and 
switches over to the stable ground stated. At T > T* 
the SC metastable phase can cross over to the stable 
ground state with a probability (p) which is governed by 
the Arrhenius equation 

p oc exp(-U/k B T) (3) 

which tells us that the barrier will be crossed with an 
ensemble average time constant r oc 1/p. If we imagine 
an ensemble of such SC phases with the same barrier U, 
then the volume of the metastable phase will exponen- 
tially decay with a time constant r. 

The energy barrier U that we talked about in the pre- 
vious paragraph is an extensive quantity. Our system is 
polycrystalline, ie., it is made up of tiny crystallites of 
different sizes. Each one of the crystallites will have an 
energy barrier proportional to its volume. Thus for a 
crystallite the energy barrier separating the metastable 
state from the ground state can be written as 

U (xVf(T ~T*) (4) 

where V is the volume of the crystallite. As already men- 
tioned / is a continuous function which vanishes for non- 
positive values of its argument. This means that various 
crystallites with the same T* will have different energy 
barriers depending on their size which implies that the 
time constant r, of the previous paragraph, will spread 
out and become a distribution of time constants depend- 
ing on the distribution of the size of the crystallites. This 
can give rise to the volume of the metallic state decaying 
in a stretched exponential manner with tim o 16 ' 17 . This 
behavior of the metallic volume with time will give rise 
to the resistivity also evolving with time in a similar fash- 
ion. We shall see later in the discussion how the metallic 
volume, the resistivity, and their time dependences are 
connected with each other. 

Imry and Wortis have argued that a weak disorder 
can cause a distribution of T*'s in different regions of 
a sample^^. This has been experimentally observed 
in a vortex lattice melting experiment in a high Tc 
superconductor^. In our case grain boundaries, shape 
of the crystallite etc. could be sources of disorder which 
would give rise to different T*'s for different crystallites. 
This effect will further broaden the distribution of time 



5 



constants in our system which arise from the distribution 
of crystallite sizes. 

Coming back to the heating runs we note that the lack 
of time dependence of the resistivity data indicates that 
the system is in a stable state. Now the question arises: 
what is the reason for the width of about 10 K observed 
in the insulator to metal transition? Could it be a super- 
heating effect? It cannot be because, as already noted, 
the lack of time dependence during heating runs rules 
out the possibility of any metastable phase in the sys- 
tem. The broadening while heating is most likely caused 
by the rounding of the phase transition due to finite size 
of the crystallites. A similar observation on broadening 
of phase transitions due to finite size of grains have been 
reported in a CMR manganite2£. 

Now we understand the reason for the hysteresis seen 
in this system. During a heating run from low tempera- 
ture the system is insulating and hence its resistivity is 
high. While cooling some of the crystallites transform 
themselves into the stable insulating state while others 
remain in the supercooled metallic state . The presence 
of the metallic crystallites lowers the resistivity of the 
system. Thus we get different values of resistivity dur- 
ing heating and cooling runs giving rise to hysteresis in 
resistivity. 

The picture so far: NdNiC>3 is polycrystalline and thus 
consists of tiny crystallites. Below Tmi, while cooling, 
the tiny crystallites can be in either a supercooled metal- 
lic state or a stable insulating state. A crystallite in the 
supercooled metallic state can make a transition to the 
stable insulating state if it crosses the barrier U with the 
help of energy fluctuations which in our case are ther- 
mal in origin. Each crystallite would undergo the tran- 
sition independent of each other. This is supported by 
the fact that the resistivity, and as we shall see later, 
the volume of the metastable metallic state, behave in a 
stretched exponential manner with time. At a particular 
temperature, as a function of time, the transition from 
the metallic to the insulating phase will begin with the 
smaller crystallites, because the barrier is small for them, 
and then proceed with the larger crystallites, in the in- 
creasing order of barrier size. While heating the system 
remains insulating and stable all the way up to the M- 
I transition and then it goes over to the stable metallic 
state. There is no superheating during heating runs. 

Now that we have understood the origin of time de- 
pendence and hysteresis in the system we would like to 
consider some more issues that needs to be understood. 
These are 

1. The dependence of resistivity on cooling rate as 
shown in the inset of Figure [TJ 

2. The behavior of the minor loops seen in Figure [TJ 

3. The intermediate ageing behavior seen in Figure [3j 

4. Both pi/po and r of Equation ([2j) go through a peak 
around 145 K and 147.5 K respectively as shown in 



Table [TJ and also in the inset of Figure [3J and Figure 

El 

Let us take up these issues one by one. 

A. Cooling Rate Dependence 

In the inset of Figure [TJ we saw that while cooling a 
lower cooling rate increases the resistivity of the sample 
in the hysteresis region. A lower cooling rate allows a 
larger number of the metastable metallic crystallites to 
switch over to the insulating state and thus the resistivity 
will be larger if the measurement is done slowly. 

If the measurement is done infinitely slowly in the cool- 
ing run we should get the uppermost curve shown in the 
inset of Figure [TJ It is interesting to note that this curve 
does not come anywhere near the heating run curve. To 
paraphrase, even if we sit for an infinite time on the cool- 
ing curve we cannot reach the heating curve. The reason 
for this is that only a small number supercooled crystal- 
lites which have a relatively small energy barrier are able 
to switch over to the insulating state at any given temper- 
ature. The others with the larger energy barriers remain 
trapped in the metallic state. These crystallites with the 
larger barriers can switch their state only if the temper- 
ature is lowered such that T comes sufficiently close to 
T* and thus the barrier becomes small and easy to cross 
(Equation ((U)). 

B. Minor Loops 

The minor loops a, 6, and c shown in Figure QJ follow 
the full hysteresis loop in the cooling cycle till we start 
the heating. While heating, curve a, with the lowest 
minimum temperature (140 K), shows a falling resistiv- 
ity with increasing temperature. Curve b with the inter- 
mediate minimum temperature (146 K) shows a rather 
flat resistivity with increasing temperature, while curve 
c with the highest minimum temperature (160 K) shows 
a rising resistivity with increasing temperature initially 
before falling and joining the full hysteresis loop at the 
MI transition temperature 195 K. 

The behavior of the minor loops can be understood 
by recognizing that the nature of resistivity, when heat- 
ing is resumed, in the metastable region is determined 
by the competition between the temperature dependence 
and time dependence of resistivity. In the case of curve 
a, with the lowest minimum temperature, most of the 
crystallites will be in the insulating state as can be in- 
ferred from the large value of resistivity. Only few of 
the crystallites will be in the metastable state and hence 
the time dependence effects can be expected to be small. 
Thus during heating the resistivity of the sample will be 
dominated by the temperature dependence of resistivity 
of the insulating crystallites. This is the reason the heat- 
ing part of the curve looks more or less parallel to the 
main heating curve before the onset of the transition. 
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The behavior of curve c is easier to understand and 
hence let us consider it before discussing curve b. In this 
case the resistivity is low and most of the crystallites are 
in the metallic state. The metallic crystallites form a per- 
colating network and the resistivity of the sample will be 
governed by their behavior. The insulating crystallites 
are essentially shorted out and will have little role in de- 
termining the resistivity. The metallic crystallites, being 
metastable, will switch to the insulating state governed 
by Equation J3]) , thus slowly increasing the resistivity of 
the sample as a function of time. When the resistiv- 
ity values are collected while heating, there will be some 
time difference between two data points depending on 
the heating rate, allowing a rise in resistivity to register. 
We would then find that the resistivity increases with in- 
creasing temperature. Curve 6 lies between curves a and 
c and in this case there will be a good number of both 
insulating and metallic crystallites. Here there would be 
competition between the temperature dependence of the 
insulating crystallites and the time dependence of the 
metallic crystallites and hence one can expect that the 
curve b would be more or less flat at the onset of heat- 
ing. 



C. Intermediate Ageing 

Now let us consider the intermediate ageing data 
shown in Figure [31 We note that after ageing for one 
hour, when the cooling is resumed, the slope of the the 
p vs. T curve is less in magnitude than for the curve 
obtained without ageing. When we stop the cooling and 
age the sample at a fixed temperature the supercooled 
metallic crystallites with relatively small U will switch 
over to the insulating state. This decreases the number 
of crystallites in the metallic state and raises the resis- 
tivity. So now, when the cooling is resumed, the number 
of metallic crystallites which can switch to the insulating 
state would be less and hence the change in resistivity 
with temperature would be lower. 

From Figure[3]it is clear that on resuming cooling after 
an intermediate stop of one hour the cooling curve merges 
with the curve obtained without ageing within about 3 K 
or less. As already noted in the previous paragraph when 
we stop the cooling and age the sample at a fixed tem- 
perature the supercooled metallic crystallites with rela- 
tively small U will switch over to the insulating state. As 
can be inferred from Equation Q those crystallites will 
have a small U which have (i) their metastability tem- 
perature (T*) close to their temperature (T) (ii) a small 
size. The major contributions to resistivity change will 
come from the relatively larger crystallites making the 
transition from the metallic state to the insulating state. 
In the light of this the merger of the cooling curves with 
and without ageing within a small temperature change of 
about 3 K means that most of the crystallites which un- 
dergo the transition from metal to insulator during age- 
ing have their T* within a few kelvin of the temperature 



of the sample. We had already reached this conclusion 
earlier when we discussed the cooling rate dependence. 



D. 



Peaks in pi / po and r 



We had noted earlier that both p\ / po and r go through 
peaks around 145 K and 147.5 K respectively. The peak 
in pi /pa is not surprising because it shows up in a re- 
gion of temperature in the cooling curve where the re- 
sistivity changes very fast with temperature. Thus one 
would expect that a large number of metastable crystal- 
lites are switching over to the insulating state near this 
region. Hence if a relaxation measurement is done here 
we should see the effect of a large number of crystallites 
transitioning to the insulating state as a peak in pi/po- 
But the peak in r at about 147.5 K is not that easily 
understood. To check whether this peak has anything 
to with percolation we should examine how the metallic 
and insulating volumes depend on temperature. 

The conductivity of a phase separated system depends 
on the volume fraction of insulating and metallic phases, 
their geometries, distribution, and their respective con- 
ductivities (a i and cr M ). McLachlan has proposed an 
equation based on a general effective medium (GEM) the- 
ory for calculating the effective electrical conductivity oe 
of a binary MI mixture 2 ^, which is 
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where / is the volume fraction of the metallic phases and 
A = (1 — /c)//o fc being the volume fraction of metallic 
phases at the percolation threshold, and t is a critical 
exponent which is close to 2 in three dimension a 22 i 23 . 
The constant f c depends on the lattice dimensionality, 
and for 3D its value is 0.16—. The GEM equation has 
been successfully applied to a wide variety of isotropic, 
binary, macroscopic mixtures and it works well even in 
the percolation regim o 21 ' 23 ' 25 . 

In order to calculate the volume fraction of metallic 
and insulating phases from Equation ((5]), we need their 
respective resistivities pm(T) and pi{T). pm{T) was ob- 
tained using pM = po + AT, where A is the tempera- 
ture coefficient of resistivity estimated from the resistiv- 
ity data above the M-I transition. pi(T) was calculated 
using the parameters obtained by fitting the resistivity 
data below 115 K to Equation (TJ. Using the above in- 
formation and the resistivity data, we calculated the vol- 
ume fraction of metallic and insulating phases and it is 
shown in Figure [H In the cooling cycle the volume frac- 
tion of the insulating phases ( V = 1 — /) slowly increases 
on decreasing the temperature below Tmi, while in the 
case of heating cycle it remains nearly constant up to 
about 185 K, and then drops to zero by 200 K. As can be 
seen from Figured] the percolation threshold for the cool- 
ing runs occurs at around 144 K. Below this temperature 
there will be no continuous metallic paths in the system. 
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Figure 4: (Color online) Temperature variation of V (= 1 — /), 
the insulating volume fraction. The blue circles show the cool- 
ing cycle (2 K/min cooling rate) and the red squares represent 
the heating cycle. The green pluses represent the insulating 
volume fraction for 0.2 K/min cooling rate. The dashed hor- 
izontal line at 84% insulating volume represents the percola- 
tion threshold. The inset shows the variation of \dV/dT\ (blue 
circles: cooling rate 2 K/min, and green pluses: cooling rate 
0.2 K/min) and AVoo (black squares) with temperature. 



The rate of change of metallic volume fraction 
(—dV/dT) has a maximum around 147.5 K in the cool- 
ing runs (inset of Figure H|). Figure [H displays the in- 
crement (AV) of insulating volume fraction with time, 
that has been extracted from the time dependent re- 
sistivity data of Figure [21 Just as in the case of time 
dependence of resistivity, the AV(t) curves of Figure [U 
were fitted to a stretched exponential function AV(t) = 
AVoo(l — exp(— (t/r) 7 )), where the constant AVoo gives 
the increase in the insulating volume fraction when wait- 
ing for an infinitely long time. The inset of Figure [4] 
also shows the variation of AVoo with temperature, which 
mimics the behavior of \dV/dT\. Figure [6] compares the 
constant of relaxation r obtained from the fitting of p(t) 
and AV(t) curves of Figure [2] and [5] respectively. For the 
p{t) curves, r has a peak around 147.5 K while for the 
AV(t) curves, surprisingly, r is nearly constant. 

During cooling \dV/dT\ of the inset of Figure |4] rep- 
resents the amount of volume that will change from the 
supercooled metallic to the insulating state for a unit 
temperature change. Now a small change in temperature 
will result in those supercooled crystallites with their T* 
falling in that temperature range switching to the insu- 
lating state. This means that the value of \dV/dT\ at T is 
a good measure of the fraction of supercooled metastable 
crystallites which have their T* close to T. Thus \dV/dT\ 
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Figure 5: (Color online) Increase in volume fraction of insu- 
lating phases (AV) with time, at various temperatures. 
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Figure 6: (Color online)Temperature variation of r obtained 
from the stretched exponential fitting of p(t) (red circles) and 
AV(t) (black squares) curves. 



of the cooling curve represents the distribution of T*'s 
in the system. In coming to the above conclusion we 
have disregarded the small fraction of crystallites that 
would be switching due to the time elapsed in covering 
the small temperature change. This is supported by the 
fact that the \dV/dT\ values calculated from the 2 K/min 
and 0.2 K/min cooling curves are essentially the same. 

The peak in AT4o occurs at 147.5 K which is also where 
the distribution of T*'s peak. This is understandable be- 
cause all the crystallites with T* close to 147.5 K would 
be switching to the insulating state if we wait for a suf- 
ficiently long time. Since the maximum number of crys- 
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tallites with T* are close to 147.5 K, we will get the max- 
imum temperature dependence for AV^o at 147.5 K. 

One would have expected that the peak in pi/po and 
AV^ would occur at the same temperature. The peak 
in pi/po is at 145 K which is 2.5 K below the peak of 
AV^. The reason for this is that the resistivity would 
change faster for a given metallic volume change nearer 
the percolation threshold. We guess that even if the vol- 
ume change is smaller at 145 K a larger resistivity change 
due to the proximity to the percolation threshold com- 
pensates for the smaller volume change and shifts the 
maximum of pi/po to that temperature. 

As can be seen from Figure [6] the parameter r of the 
resistivity time dependence fit has a peak around 147.5 K 
while there is no such peak in the r of the volume fraction 
time dependence. The rather flat temperature depen- 
dence of the r of the volume fraction is probably an indi- 
cation that the distribution of energy barriers look more 
or less the same at all temperatures. We note that the 
difference in behavior of the t's from resistivity and vol- 
ume fraction is most pronounced close to the percolation 
threshold, while away from the threshold both t's be- 
have in a very similar manner. In a resistivity relaxation 
measurement, close to the percolation threshold, the re- 
sistivity will change significantly even when the metallic 
volume changes very little. Thus, one will find that, near 
the percolation threshold, the resistivity will relax over a 



longer time period. 



V. CONCLUSION 

Our experimental results demonstrate that while cool- 
ing the physical state of NdNi03 is phase separated be- 
low the M-I transition temperature; the phase separated 
state consists of SC metallic and insulating crystallites. 
A metastable metallic crystallite switches from the metal- 
lic to the insulating state probabilistically depending on 
the closeness of the temperature of metastability and the 
size of the crystallite. At low temperature, below 115 K 
or so, the system is insulating, all the crystallites hav- 
ing switched over to the insulating state. While heating 
the crystallites remain in the stable insulating state till 
the MI transition temperature is reached and then they 
switch over to the metallic state. 
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